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Abstract 

J3 Within the framework of relativistic fluctuating hydrodynamics we compute the contribution 

of thermal fluctuations to the effective infrared shear viscosity of a conformal fluid, focusing on 
ch \ quadratic (in fluctuations), second order (in velocity gradients) terms in the conservation equations. 

Our approach is based on the separation of hydrodynamic fields in soft and ultrasoft sectors, in 
which the effective shear viscosity arises due to the action of the soft modes on the evolution 
of the ultrasoft ones. We find that for a strongly coupled fluid with small shear viscosity-to- 
entropy ratio n/s the contribution of thermal fluctuations to the effective shear viscosity is small 
but significant. Using realistic estimates for the strongly coupled quark-gluon plasma created in 
heavy ion collisions, we find that for rj/s close to the AdS/CFT lower bound l/(4ir) the correction 
is positive and at most amounts to 10% in the temperature range 200-300 MeV, whereas for larger 
values n/s ~ 2 /(Air) the correction is negligible. For weakly coupled theories the correction is very 
small even for n/s = 0.08 and can be neglected. 
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I. INTRODUCTION 



The heavy ion collisions experiments performed at RHIC and LHC create a hot and 
dense medium, the so-called Quark-Gluon plasma (QGP), that is believed to be strongly 
coupled. The most compelling evidence supporting this idea comes from relativistic viscous 



fluid dynamics simulations that re 



by experiment (see, e.g. Refs. jl4l2j). The momentum anisotropy is the translation to 



3roduce the momentum anisotropy patterns measured 



momentum space of the initial spatial eccentricity of non central collisions. A weakly in- 
teracting system of particles can not convert spatial anisotropy into momentum anisotropy 
in an efficient way, but this translation can occur quite efficiently if the particles interact 
strongly. The observed large magnitude of the momentum anisotropy indicates that the 
QGP is strongly coupled. Another indication of the strongly coupled nature of the QGP is 
the strong quenching of high-p^ probes measured at RHIC l|-|4j]. 

Transport coefficients such as the shear (rj) and bulk (() viscosities are crucial inputs in 
fluid dynamics simulations attempting to describe the evolution of matter created in heavy 



ion collisions. Great efforts are currently focused on developing new theoretical 
compute these transport coefficients more accurately from microscopic models 



also on extracting them more precisely from RH 



obtained from Lattice QCD calculations 
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C and LHC measurements 
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tools to 
and 



. Results 
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301 ] indicate that rj/s and (/s, where s 



is the entropy density, depend significantly on temperature, showing a minimum and maxi- 
mum, respectively, at the critical temperature corresponding to the QGP-hadron crossover. 
Note that although it is expected that the transport coefficients of the QGP depend on 
temperature, the impact of a temperature-dependent i)/s on momentum anisotropies as 



obtained from simulations has been investigated only recently 
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35|. 



It is possible to extract an average value of rj/s by matching the particle spectra and 
elliptic flow obtained from fluid dynamics [8| or hybrid fluid-kinetic simulations to 



data. See also Refs. [36|-|38[ for a description of other approaches to extract values of 
rj/s from data. The optimal value of rj/s that comes out of all these fits is close to the 
lower bound rj/s = l/(47r) that was found by Kovtun, Son and Starinets [39| (KSS bound 
from now on) via the Anti de Sitter/ Conformal Field Theory (AdS/CFT) correspondence. 
However, deviations from the KSS bound by a factor of two or more are still possible due to 
uncertainties stemming from diverse sources such as the initial conditions, the dynamics of 
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chira 
5-12 



nc 



35 



ds coupled to the quark fluid, and the freeze out stage, to mention a few - see Refs. 



40 



45| . Very low values of rj/s are also favored by studies of the overall entropy 



production during a collision - see e.g. 46] - and, albeit more indirectly, by the fact that 
flow structures such as cones and ridges are able to survive until freeze out. If dissipative 
effects were too large, such collective flow structures would be washed away much earlier in 
the evolution of the fireball 
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50|. 



The system created in heavy ion collisions consists of a fireball of quarks and gluons that 
expands and cools very rapidly under its own pressure. Since the fireball is formed out of a 
finite number of particles and has a small size in the range 5-10 fm, it is natural to expect 
that thermal fluctuations may have an observable impact on its evolution. In spite of this 
expectation, the role of thermal fluctuations and its impact on hydrodynamic evolution has 
only very recently been discussed in the context of heavy ion collisions 51n53 ] . 



It was shown in Ref. [52] that for the QGP with small values of t]/s near the KSS bound, 
thermal fluctuations on top of long-lived sound and shear waves can contribute significantly 
to the effective value of 7]/s as would be measured on large scales. In contrast, for larger 
values rj/s ~ 2/(47r) the correction coming from thermal fluctuations on sound and shear 



waves was found to be negligible. Specifically, in [52] the authors used Kubo's formula 
together with second order hydrodynamics and focused on corrections to f]/s coming from 
the zeroth order (in gradients), quadratic term (p + Pq)u^u u in the energy momentum 
tensor T M1/ , where po and Pq are the equilibrium energy density and pressure (Pq = po/3 
for a conformal fluid), respectively, and is the flow velocity fluctuation ( 52] considers 
a vanishing background velocity). A related study on the effective shear viscosity of the 
strongly coupled QGP was carried out in [53] using the AdS/CFT correspondence, where 
corrections to r] arising from higher order velocity gradients were computed using the idea 
of resumming those corrections into an effective t](u, k) depending on both frequency and 
momentum. 

In this work we compute the correction to the shear viscosity that come from the effect 
of relatively small scale fluctuations of the hydrodynamic modes on their evolution on larger 
scales. For this purpose, we use the equations of relativistic fluctuating fluid dynamics 
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58] and focus on terms which are second order in velocity gradients (in T^ v these 



terms are first order) and quadratic in velocity fluctuations. The development presented here 
is based on a division of hydrodynamic modes in ultrasoft and soft sectors, both of them 
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subjected to the influence of stochastic noise coming from hard (i.e. particle-like ) m odes. It 
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and 



is, therefore, inspired in the hard thermal loops approach to thermal field theory 
particularly in earlier studies of the Langevin dynamics of soft and ultrasoft modes 
As this approach does not rely on Kubo's formula it is different from the one adopted in 
We anticipate that the correction computed here is of order \ 2 = PmaxV 2 / ( S T) 2 or higher, 
where p ma x is the value of momentum beyond which the second order gradient expansion 



67]. 
52]. 



breaks down, while the correction obtained in [52| is of order \. Since \ must be small 
for a fluid description to be reliable, the correction to rj/s obtained here is suppressed with 
respect to the one of 52|, but as we will show it is still significant. We will compare our 
results with those of [52] in Sections II III and IIII Al 

It is appropriate at this point to emphasize that our approach is purely phenomenological. 
It is not our intention to give a precise value of the correction to the microscopic shear 
viscosity arising from second order terms, but rather to provide a reasonable upper bound 
for it under realistic conditions for the QGP. We note that the correction that we obtain, 
as well as that obtained in 52|, is inversely proportional to rj/s, and therefore diverges for 
vanishing rj/s. As pointed out recently by Torrieri [68| in the context of ideal quantum 
hydrodynamics, the cure for this divergence may lie in quantum corrections that set in for 
small rj/s and therefore "correct the correction". Although this is an important issue with 
strong implications for the study of matter created in heavy ion collisions and surely deserves 
further attention, here we shall not dwelve further into it since we focus on the correction 
to r]/s arising at finite values of this ratio. 

The paper is organized as follows. In Section [IT] we describe the theoretical setup and 
obtain the equations for the ultrasoft and soft modes. In Section [III] we compute the correc- 
tion to the shear viscosity arising from second order soft terms, and then provide numerical 
estimates for a weakly and a strongly coupled theory, in particular focusing on the QGP 
created in heavy ion collisions. Finally, we conclude with a brief summary and outlook in 
Section [IV] Details on the calculation of the correction to rj/s are given in Appendix lAl 
while a description of the procedure used to compute the frequency that separates ultrasoft 
and soft modes is given in Appendix [B] 
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II. FLUCTUATING CONFORMAL HYDRODYNAMICS 



A. Soft and ultrasoft hydrodynamic modes 

We consider a conformal fluid in flat spacetime with signature (+, — , — , — ). The comoving 
time and space derivatives are D = u^d^ and V M = A^c^, respectively. The hydrodynamic 
equations at second order in velocity gradients read 

ipDu"- \vp + 7iU&po>* = f» 

3 4 (1) 

Dp + -pVX - rja^a^ = f 

where we have made the approximation II M1/ = r]<j fJ-l/ = i]V <fl u u> which is enough for our 
purposes here. The brackets <> around indices denote the symmetric and traceless pro- 
jection orthogonal to w M , while round brackets denote symmetrization. 77 is the microscopic 
shear viscosity as derived e.g. by Kubo's formula, and (/, / M ) are stochastic noises. For 
recent reviews on relativistic viscous fluid dynamics as applied to heavy ion collisions see 
Refs. (5- T] . 

Next we divide the hydrodynamic fields into "ultrasoft" and "soft" parts, 

p = p< + p> and t/ 1 = u< + m> (2) 

On top of this separation into ultrasoft and soft sectors, we perfom an expansion in deviations 
from equilibrium characterized by a small parameter, e, staying at linear order, thus 



P< = Po,< + epi,< + • • • 
P> = Po,> + epi,> + • • • 



u < = u o,< + eu i< 



u > = u o,> + eu i,> + • • • 



(3) 



In this way (po,< + Po,>;^o< + u o>) correspond to thermal equilibrium. For notational 
simplicity is it convenient to define 

Po = Po,< , Pi = Pi,< 
$Po = Po,> , Spx = pi,> 

(4) 

V* = <,< , = Ulj< 
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It will be shown in the next section that V 1 = (V°, 0, 0, 0), with V° close to but not exactly 
unity, and that = V^W^ = 0, which imply that s M = (0, s) and W = (0, W) in the 
local rest frame (LRF). We note that ^ 0. Additional constrains arising from the 
normalization of u^ 1 are discussed in the next section. 

Note that our treatment is purely hydrodynamic, i.e. assuming momenta to be small, so 
no attempt is made to actually derive the noise as coming from the hard sector of a micro- 
scopic field which should be described by transport equations such as Boltzmann equation 



56l . I63I-I67]. So, the ultrasoft and soft modes of the hydrodynamic fields are described by 



the conservation equations of fluid dynamics, while the effect of the hard modes on the 
hydrodynamic evolution is included through the noises (/, This is not a limitation to 
our purpose here, since we intend to determine the correction to rj coming from long-lived 
hydrodynamic modes (shear and sound waves). 

A reasonable estimate for the value of the momentum separating soft and hard modes will 
be given in Section IIHI This estimate is essentially based on the requirement that second 
order terms in derivatives must be smaller than first order terms, i.e. that the gradient 
expansion does not break down. The value of momentum that separates ultrasoft from soft 
modes can^be computed by requiring that the equations of the soft modes be linearizable 
(see, e.g., 



56 



591]). This is done in Appendix [B] Due to the phenomenological nature of 



our approach, and taking into account that we intend to provide an upper bound to the 
correction to r]/s, we also compute this correction for different values of the soft-ultrasoft 
separation frequency. We will come back to this point in Section IIHI For the moment, let 
us leave this separation implicit. 

In what follows we will use Latin indices k) to denote spatial coordinates, i.e. s M = 
(0, s l ) with % = 1,2,3. At zeroth order in e the conservation equations for the ultrasoft 
modes read 

(SD Spo) + ~ (5 Po V m s l ) - \ (paVaSi&s*) - V = 

where (•) is a thermal average. The zeroth order linearized equations for the soft modes are 

4 

D 5p + -p V( )iS J = / 

4 1 ■ k ■ (6) 

-p D s* - -V\ 0) 5p + v^\o)A sJ = f ■ 
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Note that the noises are determined by the equilibrium state, as dictated by the fluctuation- 
dissipation theorem. Their role is to enforce thermal equilibrium at zeroth order in e, 
constraining the expression for the thermal correlator of velocity fluctuations which will be 
given in Section ITHl 

At first order we get for the ultrasoft modes 

(5D 5 Pl ) + (SD 1 Sp ) + ^PoV ( o)^ + ^ (SpiV^s') 

- I (poV iflS4>) d^) + ± <5poV (0 )^) - | (8 Po V s i d°W i ) (7) 

- \ (5p V W t d°s l ) - (5p W {lSl) d l si) - ^ ( Po V t t d s*) - 2 V ( SlJ t ij ) = 



and 



+ ^ (5p D t») + ~ (5po5D W») + ~ (8p D lS ») 
+ ± (6 Pl D s») + jj (V^d^px) + 2 - (W^s 
+ l(yH^d 4> Sp )+r ] C^O 



with 



- 2g ia (wWdf,s afi ) - 2g ia (vWdps *) . 
For the soft modes we get 

D 5 Pl + ZVp - ^Po WidV 

+ ~PoV( 0)/ ^ + ^SpoV^W* - 2 VSij W ij = 



and 



with 



(8) 



(9) 



(10) 



ipoA)^ + ^p ZW M - ^Vf 0) 5 Pl + j^W^<VPo + tj£" = (11) 



E M = Af 0)Q ^t Q/3 - 2^ 7a \/ (/i s 7) ^iy a/3 - 2g 1Q V^W l) d fi s^ . (12) 
In deriving Eqs. (]5])-( |l2"]) we have used that 

V" = Vf 0) - 2V^s (t>) d <t> - 2eV^W* ) d 4> - 2eW { »s' t,) d <t> - 2eV { H 4>) 8 4> (13) 
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where Vf 0) = (g^ - VV v )d v and = (A» u + A v »)/2, 

D = D + 5D + e£>! + e5D 1 (14) 
where L> = Wdp, SD = s^, D x = Wd^ and 5D 1 = t"d^ and put 

S — V (Q) i> 

t"" = Vg^ (15) 

which by definition are traceless and orthogonal to V^. Note that we have neglected terms 
0(5 3 ) or higher in the thermal averages. Moreover, we have set p± = 0, which is a reasonable 
approximation since, as it will be seen later, the term we are interested in is already linear in 
e. Note that, since we stay at first order and 0(5 3 ), those terms containing V <fl that would 
appear in Eqs. fl5|)- f|T2|) actually vanish because in these terms we can put A^ u = (see 
Eq. (HI). 



B. Constraints on four velocity 

We will now discuss the normalization of the fluid velocity. For consistency, the velocity 
should be normalized both in the mean and at linear order in fluctuations, i.e. u^u^ = 1, 
because it is who satisfies the conservation equations ([T]). We also require that (V 1 + 
s M ) 2 = 1 which implies 

V"% + = 1 and 

(16) 



so at first order we get 



K% = 



V»W lt = and 



V% + W% = 



(17) 



The transversality and tracelessness of a^ v must be satisfied both in the mean and at 
linear order in fluctuations. The constraints coming from the transversality and tracelessness 



S 



of then read 



(s*V) = 

(yb 8 +% 8ll ) = o 

= 

+ W"V = 
V^s^d^ + V^W^d^ = 



These relations were used in deriving Eqs. (l5l)- (jT2l) and will be used in what follows. 

III. SHEAR VISCOSITY INDUCED BY THE SOFT MODES 

The idea is now to solve the soft mode equations to compute the induced viscosity arising 
in the ultrasoft equations. We will solve the soft mode equations in Fourier space. For any 
given quantity R{x^) we have 

R{W) = ^[R{x)){k") = J R{x")e ik ^ d A x (19) 

where ^[g(x)](k fJi ) = g{k^) denotes the Fourier transform. We split the wavevector = 
un^ + where = V^/\V 2 \. Since V 1 = (V°,0) in the LRF and V^s^ = we have 
s M = (0, s), and similarly for (but not for t^). Note that p^ = (0,p l ) in the LRF. 
The zeroth order equations then read 



4 

uopo + ^PoPjS 3 = if and 



4 ~i 1 i x ~ ir l 
-p us - -p 6p - — 



1 . 



if, 2 = 1,2,3 



(20) 



Without loss of generality we can set p 2 = P3 = 0. The solution is (jp = px] 

irjfp 2 + 2fp u + 2f 1 p p 



5po = 3i 
§i = 3i 



A 

3/iW + fp 



32 B 
- 6/ 3 
53 = ~B 



2A 

6/2 



(21) 



where 



A = 3ir\ujp 2 + 6po(u 2 — -p 2 ) and 



B = 3r]p — 8ipou 
The first order soft mode equations are 



(22) 



oo5p 1 + -p oPj t j = -iG[W°] 



-pXp^) + typj 



where 



d x 



2r 1 g ia V^s l) dBW a l 3 - 2r]g ja V^W^dBS 



a/3 



d 4 x 



whose solution is 



Spi = 3i 



irjGp 2 + 2Gp u + 2H 1 p p 



A 



ti = 3i 



.3HiU + Gp 



2A 



QH 2 
B 

B 



tn = %- 



Ap UJ 



(23) 



(24) 



(25) 



We are interested in those terms appearing in Eq. (fl~T|) containing two (orthogonal to 
V 1 ) derivatives of W 1 '. Moreover, since the bare term A^ a dpW a ^ is orthogonal to V^, the 
induced terms that provide the correction to rj must be orthogonal to as well. The only 
term in proportional to dpW a ^ is 



(26) 
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since V a dpW a ^ = 0. Therefore, this term in is proportional to V^. So, those terms 
in Eq. ( II ip which upon replacement of f 1 by its expression given in Eq. (125]) contain two 
derivatives of are not orthogonal to V^, but proportional to it instead. Examples of 
such term are (45poD t^ /3) and —2r]g^ a (V^P'dpS^}. Note that the latter term vanishes 
identically because V a dps al3 = 0. On the other hand, a term such as (45piD s^) does not 
contribute to the correction to 77, since 5p\ depends on Hi, whose term proportional to 
dpW a ^ vanishes (it is proportional to V\ — 0). 

Therefore, the only term satisfying all the requirements mentioned above is {Ap^SD^ /3) . 
This term is orthogonal to and is second order in spatial derivatives of W 1 *. The correction 
to T], which we denote as r) c , comes from the coefficient multiplying dpW al3 in this term. We 
focus on the p = 1 component, i.e. we only consider the x component of the momentum 
conservation equation. To be definite, we take a = 2 and /3 = 1, which means that we 
are considering the correction to 77 that comes from the term d x Vj ^W x . Note that from 
SDi = t a D a appearing in the term (4po^-Dis A '/3) we only need t , because this component 
is the only one containing dpW a ^ - see Eqs. ( 1241) and ( 1251) . Taking into account these 
considerations, the term we are interested in becomes (for simplicity we use x instead of x M ) 

A(x)^^(t°(x)D (x)s 1 (x)) (27) 

together with the equation 

P{x)t°(x) = rjs 2 (x)dW(x) (28) 

where we have defined the operator P(x) = ApqDq/3, and set dW = d x W yx . For simplicity, 
we have approximated V° = 1 as well, which implies neglecting (siS 1 ) in the normalization 
of see Eq. (fTO]) . The term given in Eq. (I27p . from which the correction to the shear 
viscosity r\ c arises, is quadratic in fluctuations and second order in spatial gradients (because 
t° is proportional to dW which is second order). To obtain 7] c we must solve Eq. ([28]) and 
plug it into Eq. (12?]) . Details of this procedure are given in Appendix |A] The expression for 
r] c that results is given in Eqs. ( lAl^ -f lATS]) and Eq. (1A24]) . 

In order to give concrete numerical estimates for r] c /s in the next section, it will prove 
convenient to express 1] C in terms of dimensionless quantities. The relevant parameters 
appearing in the expression for i] c that is given in Eqs. (1A14[) - (1A18|) and Eq. (IA24I) are 
(^max,Pmax) - the maximum momentum for the soft modes, beyond which the gradient 
expansion becomes invalid -, oj m in (the frequency separating ultrasoft and soft modes), and 
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7 = t]/(sT), where s is the entropy density. We therefore define the two dimensionless 
quantities 



X = IPr, 



(29) 



and 



UJ r , 



(30) 

where t„ is the second 



^max 

The value of p max can be roughly estimated as p ma x ~ 1/ (2r n 
order coefficient representing the relaxation time of the full shear tensor IP 1 ' towards its 
Navier- Stokes value IP" = i]a fiu . The value of r v has been computed in strongly coupled 



M = 4 super Yang-Mills theory [25 
coupled theory from kinetic theory 



-|27l| . resulting in r n ~ 2.Qi]/(sT) = 2.67, and in a weakly 



24j, resulting in 77 ~ 57. The trend is that 17/7 increases 



as the coupling decreases. Here, we shall use as a reasonable estimate a value of = 37 to 
model the strongly coupled QGP, and a value of = 67 to model a weakly coupled theory. 
It is important to note that, as it has been recently recognized 28(, the true relaxation time 
appearing in the equation of motion of dissipative currents is given by the first pole of the 
retarded Green function. It turns out that the location of this pole can not be found, in 
general (and particularly for a strongly coupled theory), from a gradient expansion. In other 
words, the equation of motion for the dissipative current is not a relaxation-type equation 
for details. Nevertheless, in this paper we will use r n oc 7 as computed from the 



sec 



second order gradient expansion 24j-|27|, both for weakly and strongly coupled settings. 

We note that by taking the value p m ax = 1/(2^), the validity of the gradient expansion, 
i.e. the fact that higher order terms are smaller than lower order ones, is ensured, since for 
this to happen one must have p ma x < I/tv (of course, second order fluid dynamics becomes 
more and more reliable the smaller p m axT-n is) - see 52) . This condition also ensures that 
hydrodynamic waves are actually long-lived excitations of the fluid and thus well defined. 
As an estimative value for 0J max , we shall take 0J max ~ l/(2r 7r ) as well. In Section [III Al we 
will consider other values of p ma x and analyze the dependence of our results on these values. 

It is useful to note that since 77 ~ rj/(sT), a term containing p max T in rj c is proportional 
to T 3 (r]/s)~ 2 . If we put s = aT 3 (estimates of the value of a for strongly and weakly coupled 
theories are given in the next section) we get that 

,2 



Pn 



SX 



a(rj/s) 



(31) 
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To simplify the resulting expression for r] c /s, let us define the following quantities 



1 + v^ + X 



M 



1 - a/2x + X 

x(i + 
i-x 2 

x " v and 



and 



X 



i-^Wx + x 



i-x 



(32) 



r 



27r 3 a(?7/s) 2 
We then have (see Eqs. (jXl4jl - (jXl8jl ) 



Vc 

s 



+ 



K 

2 ^3/2 



2X 



1 



4 1 

-(1 — 7ty — arctan M — 

3 V 3 3 3v^x 



In L r 



In L r 



arctan N„ 



4^3/2 

+ - — ;= arctan N„ 



(33) 



3^2x 3^ ""~ 3^2x 

The expression for the correction to the shear viscosity given in Eq. fl33|) is our main 
result. We have therefore computed the correction to rj that comes from the action of the 
soft modes on the evolution of the ultrasoft modes, as is reflected in the second order term 
d x V v ,$W x of the fluid dynamic equations. 

At this point, it is important to stress that out treatment is, although similar in spirit, 
otherwise quite different from the one put forward in Ref. 52| in two main points. First, 
the developments of 52) are based on Kubo relations for rj and t v in terms of the retarded 
correlation function for two T xy stress tensor operators, while ours rely on the explicit 
separation of hydrodynamic modes into ultrasoft and softs components that is done on top 
of an expansion in deviations from thermal equilibrium. Second, the authors of 52] focus on 
the term (po + Pq)u^u u of the stress energy tensor (the background velocity is zero), which 
does not contain gradients, while we focus on an induced term in the conservation equations 
that is second order in (background) velocity gradients (therefore the term comes from a first 
order term in the stress energy tensor). The contribution to r] c of the term that we consider 
is suppressed with respect to the zeroth order term by extra powers of x — Vmaxl- The 
zeroth order term ~ u 2 leads to a correction to r]/s that is of 0(x), which in our notation 
reads [see Eq. (4.2) of [521]: 

UX (34) 



Vc 

s 



I20ir 2 a(ri/s) 2 

As it can be seen from Eq. (133]) . the contribution to r) c /s found here is 0(x 2 ) or higher. 
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In the next section we will analyze the behavior of r] c /s with rj/s for strongly and weakly 
coupled theories and compare the obtained results with those of [52]. 



A. Numerical estimates for weakly and strongly coupled theories 

We will now provide some estimates for the strongly coupled QGP as well as for an 
illustrative weakly coupled theory. We shall take the entropy density to be s ~ aT 3 , with 
a = N% (N c is the number of colors) for a weakly coupled th eory and a = 10 as a typical 



value for the QGP in the temperature range 200-300 MeV (see |69|, |70|). With respect to the 
second order transport coefficient t w , for the QGP we will take tv ~ 37 [25], which implies 
that x — Pmaxl = 1/6 (if V max = l/{2r n )), while as a typical value for a weak coupling 



theory we will take \ = 1/12 24] . Recall also that we take ca 



Pn 



In the following we will focus on values of rj/s in the range of those values that can 
be extracted by matching fluid dynamic simulations to elliptic flow data - see js-isj]. In 
particular, although we will take the values of rj/s to be in the range 0.06-0.16, we will pay 
special attention to the comparison between the cases rj/s = 0.08 and rj/s = 0.16. The 
reason for this choice is that 11 Is = 0.08 is close to the KSS lower bound obtained from the 

n 

AdS/CFT correspondence [39], while rj/s = 0.16 is close to the maximum value of rj/s that 
is favored by fitting fluid dynamic simulations to elliptic flow data 

To estimate rj c / s we still need to provide values of co min . A natural request to our approach 
is that oJmin must be such that the equations for the soft modes can actually be linearized. 
This condition is analyzed in detail in Appendix [Bj Here we will just quote the values 
obtained from this analysis. For a strongly coupled theory with % = 1/6 we get from Eq. 
flBTOj) that f = u min /u max = 0.88 for rj/s = 0.08 and f = 0.72 for rj/s = 0.16. For a weakly 
coupled theory with x — 1/12 we get £ = 0.78 for rj/s = 0.08 and £ = 0.56 for rj/s = 0.16. 



With the estimates given above we can compute rj c /s as a function of rj/s in the range 
of values of relevance to heavy ion collisions. Figure [1] shows the percentual correction 
{rjc/s) x 100 /{rj/s) as a function of rj/s, for \ = p m ax1 = 1/6 (corresponding to the strongly 
coupled regime). The results are shown for two values of £ = uj m i n /uj max = 0.88 and £ = 0.72 
corresponding to the values of u min computed in Appendix [Bl for rj/s = 0.08 and rj/s = 0.16 
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FIG. 1: (Color online) Percentual correction (rj c /s) x 100/ (r]/s) as a function of r\js for £ = 
Wmin/oJmax = 0.88 and £ = 0.72. The value of x = Vmaxl is 1/6, corresponding to a strongly- 
coupled theory. 

(see Eq. (IBlOj) ). It is seen that for 77/s ~ 0.08 the correction amounts to roughly 8% of 
f]/s. Due to the inverse proportionality of rj c /s on the square of the microscopic ratio, the 
correction drops to zero quite fast. For rj/c = 0.16, the correction is less than 1%. 



Going over to a weakly coupled theory, Figure [2] shows, as an illustrative case, the per- 
centual correction (r] c /s) x 100/(?7/s) as a function of rj/s for x — Vmaxl = 1/12 and a = 9 
(corresponding to N c = 3 colors). It is seen that the correction r] c /s is much smaller than 
the one obtained in a strongly coupled setting. For rj/s = 0.08, the correction is ~ 2%, 
whereas for 77/s = 0.16 it amounts to less than half percent. 

The results presented so far are based on an estimation of p ma x, the maximum value of 
the momentum for which the gradient expansion is valid, as l/(2T n ). Since this estimation, 
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FIG. 2: (Color online) Percentual correction (rj c /s) x 100/(rj/s) as a function of r\js for £ = 
Wmin/oJmax = 0.78 and £ = 0.56. The value of x = Vmaxl is 1/12, corresponding to a weakly 
coupled theory. 

although reasonable, is not rigurously justified, it is important and necessary to quantify 
the impact of different choices for p max on the values of r\ c j s that we obtain. To assess the 
strength of the conclusions that can be extracted from our results with respect to the specific 
value of Pmax, Figure [3] shows r) c /s as a function of 77/s for three different values of p ma x = 
1/(717"^), namely n = 4, n = 2 and n = 4/3, with £ = co m i n /u max = 0.88 (corresponding 
to 77/s = 0.08 as computed in Appendix IB!) and x — Pmaxl = 1/6 (corresponding to the 
strongly coupled case). It is important to recall that the second order gradient expansion 
ceases to be valid for p max larger than l/r n . From the results shown in Figure [3] it is seen 
that the correction to rj/s increases with increasing values of p m ax, but the important point 
to note is that even for a value of p ma x quite close to l/r n (the case with n = 4/3) the 
correction at 77/s ~ 0.08 remains ~ 10%. Figure [3] shows that, for realistic estimates of the 
parameters corresponding to the strongly coupled QGP in the range T = 200-300 MeV, the 
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FIG. 3: (Color online) Percentual correction (rj c /s) x 100/(r//s) as a function of r]/s, for three 
different values of p max = l/(nT n ), with n = 4, n = 2 and n = 4/3. The results correspond to 
£ = Wmin/wmm = 0.88 and x = PmaxJ = 1/6 (strongly coupled theory). 

correction to rj/s coming from second order terms in the fluid dynamic equations is at most 
~ 10%. 



As discussed in Appendix El here we estimate the values of u min for different values of rj/s 
by requiring that the equations describing the evolution of the soft modes can be linearized. 
In view of the fact that our approach is purely phenomenological, it is also important to 
quantify the effect that changing the value of £ has on the results that we obtained. To this 
end, Figure H] shows r] c / s as a function of r\js for £ = oj m i n /u ma x = 0.88, corresponding to 
the estimate for rj/s — 0.08 obtained in Appendix [Bj together with two other close values 
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FIG. 4: (Color online) Percentual correction (rj c /s) x 100/(r]/s) as a function of r\js for £ = 
Wmin/oJmax = 0.88 (corresponding to the estimate for rj/s = 0.08 obtained in Appendix [B]), and 
two other close values £ = 0.65 and £ = 0.95. The value of \ = p ma xl is 1/6, corresponding to a 
strongly coupled theory. 

£ = 0.65 and £ = 0.95 chosen to illustrate the behavior of r] c with £. The results correspond 
to the strongly coupled case with \ — Pmaxl = 1/6. It is seen that the value of r\ c does not 
depend strongly on £; the variation of r] c /s with £ is small for rj/s = 0.08 and negligible for 
rj/s = 0.16. These results show that the correction to r]/s for the strongly coupled theory 
remains bounded, giving strong support to the conclusion that r] c /s < 10% for the QGP 
under realistic conditions in the temperature range T = 200-300 MeV. 

Our results clearly illustrate the fact, already highlighted in 52J, that for a strongly 
coupled fluid such as the QGP formed in heavy ion collisions, thermal fluctuations actually 
do change the value of the effective infrared r]/s, provided that the microscopic value of the 
ratio is small. Specifically, we have found that in strongly coupled theories with an 7]/s close 
to the KSS value the correction due to second order terms in the fluid dynamic equations is 
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at most ~ 10% of the "bare" value, while for larger values of rj/s it is extremely small and 
can be safely neglected. For a weakly coupled theory the correction turns out to be much 
smaller than at strong coupling, and can be neglected even at rj/s = 0.08. These are the 
main results of this work. 

To put our findings into perspective, it is appropriate at this point to note that the 



correction r] c /s found in Ref. 52J from the zeroth order term (po + Po)u^u 1 ' in T^ v is, for the 
same values of p max , s, etc. as used here, roughly 50% for rj/s = 0.08 (for rj/s = 0.16 it is 
negligible). The correction computed here, that comes from the second order ultrasoft term 
in the fluid dynamic equations, is at most five times smaller for rj/s = 0.08. Our results 
provide further support to the conclusion that, for the strongly coupled QGP created in 
heavy ion collisions, a value of rj/s = 0.08 at T = 200-300 MeV seems to be ruled out by 
the effect of thermal fluctuations on the hydrodynamic evolution of the longest wavelength 
modes, but a value of rj/s = 0.16 is still quite plausible. 



IV. CONCLUSION 



In this work we compute the correction to the microscopic shear viscosity (i.e. the one 
computed from Kubo formula) induced from thermal fluctuations appearing in second order 
gradient terms in the equations of relativistic fluid dynamics. Our calculations are based 
on the idea of separating the hydrodynamic variables in ultrasoft and soft sectors, and 
determining the effect of the latter on the evolution of the former. The hard modes, which 
are not described by fluid dynamics, enter the conservation equations through stochastic 
noises. In this sense, the developments presented here fall naturally into the relativistic 
fluctuating hydrodynamics framework. 

We find that the correction to rj/s is positive and inversely proportional to the square of 
this ratio. For the conditions prevailing in heavy ion collisions, we estimate the maximum 
correction to the viscosity-to-entropy ratio of the QGP at T = 200-300 MeV to be ~ 10% 
for rj/s = 0.08 and ~ 1% for rj/s = 0.16. These results indicate that for small rj/s (close 
to the KSS bound) the influence of soft modes on the evolution of the ultrasoft modes is 
significant. For larger values of rj/s, the corrections are indeed very small and can be safely 
neglected. The correction for a weakly coupled theory is much smaller than the one found 
at strong coupling and can be neglected even for value of rj/s close to the KSS bound. 
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Our results support the conclusion found in 52] that thermal fluctuations have a sig- 
nificant impact on the value of the effective infrared rj/s corresponding to the quark-gluon 
plasma created in heavy ion collisions, provided that the microscopic value of rj/s is small. 



Our work, and those of 



51 



lH53| . evidence the important fact that what could actually be 
extracted from data (for example by matching computed hadron observables to measure- 
ments) is an effective value of the ratio rj/s rather than the value computed microscopically. 
In this sense, our work gives further support to the conclusion that a value of rj/s = 0.08 
for the strongly coupled QGP that is created in heavy ion collisions seems impossible due to 
the effect of thermal fluctuations on the hydrodynamic evolution in the range T = 200-300 
MeV, whereas a value of rj/s = 0.16 is not ruled out. 

Probably one of the most interesting things to try out is to simulate the three-dimensional 
evolution of the fireball created in heavy ion collisions using the equations of fluctuating 
hydrodynamics in order to be able to extract a value of rj/s which in principle should be 
closer to the real one (first steps in this direction have been taken in |5JJ with a boost- 
invariant Bjorken model). Moreover, a comparison with the values extracted from viscous 
hydrodynamic simulations without thermal fluctuations would allow us to directly quantify 
the correction to rj/s arising from the effects of such fluctuations on the hydrodynamic 
evolution of matter created in heavy ion collisions. 



Appendix A: Details on the calculation of the effective shear viscosity 

In this appendix we describe in detail the calculation of the correction r\ c to the shear 
viscosity. We start by solving Eq. (128]) . thus obtaining 

t°(x) = J dy[P}-\x - y)s 2 (y)dW(y) (Al) 

so that 

A (x) = ^ f dy[P]-\x - y)D (x)S 12 (x - y)dW{y) (A2) 



3 

with Si 2 (x — y) = (si(x)s 2 (y)). Recall that P(x) = 4poD /3. The correction to r\ comes 
from the term multiplying dW, that is to say, we would like to write A(x) as r] c (x)dW(x) 
for some r] c (x). To achieve this and proceed further, we will assume that dW(y) is a very 
slowly varying function, which allows us to approximate A(x) ~ i] c (x)dW(x) with 

4p ?7 



Vc{X) 



J dy[P]-\x - y)D (x)S 12 (x - y) (A3) 
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The approximation that dW be a slowly-varying function is appropriate for our aim of 
providing an estimative upper bound on r/ c , and in any case it can be regarded as a definition 
of 77 c that is valid for all practical purposes. 
Going over to momentum space, we can write 

VeW) = j d*\ -j^S^ - A") (A4) 



or, putting = — A M , 



It is important to remark that the frequency integrations in Eqs. ()A4j) and AA5D run over 
all (positive and negative) values of the frequency Since Su is even in a (it must due to 
the invariance of the equilibrium state), we can change the integration limits in Eq. f ]A5j) 
to o~° G (0,oo). An immediate consequence of this is that the terms oc k° drop out when 
integrated, and hence need not be considered further in what follows. 

Assuming that the noises satisfy the fluctuation-dissipation theorem, the correlator of 



velocity fluctuations can be readily computed 52J, 15414581] . It is given by 



121 >~ (po + J Po)K) 2 + 7V 4 ^3(p + Po)[( C T ) 2 -^/3] 2 + (47^ /3) 2 1 ' 
with cr = o t (j i and 7 = 77/ (p + Pq) = f]j sT (s is the entropy density). The first and second 
terms correspond to shear and sound waves, respectively. 

We will compute the contributions to i] c coming from shear and sounds waves separately. 
Let us start with the shear waves. We get 

where we have defined b = (7/a ) 2 . The rationale for introducing in Eq. (1A7I) the upper 
and lower bounds on the frequency a is discussed below. 

As mentioned in Section [Til we assume that only soft and ultrasoft momenta contribute 
to r/ c , since these are the modes actually described by fluid dynamics. This means that the 
integral must be done over modes satisfying {a , a 1 ) < {u max ,p max ), where {u max ,p max ) is 
the value of the four momentum separating soft and hard modes as discussed in Section [TTJ 
The reason for introducing co m i n 7^ in Eq. ( 1A7I) is to ensure that the equations for the soft 
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modes can be linearized, as discussed in detail in Appendix |B] Estimative values of both 
cutoffs Umax and uJmin are discussed in Section IIII Al 

We now return to Eq. (IA7|) . Assuming momentum isotropy the spatial integral becomes 



/ 1 r^ dcT = o 7=1 — Id t + 2arctan 



We shall obtain the general expression for r\ c for k ^ 0. Defining 



# 



2 3/2 7r 3 7 l/2 

the first contribution becomes 



In terms of $ = a/p max we have 



and 



„(3) - 

/c 2 5 / 2 



Performing the integration we get 



3 



(2) = T lPmax 



1 _|_ $4 I max 1 i — ./9$ _J_ $2 

1 mm / \ 1 V ^^max ' ^ max 



>2 

■ max 



2 • ; 
mm J 



(A8) 



with 

^1,2 = pL ± V2ap maa; + a 2 (A9) 

and 



a = yi/b = (A10) 
o 



7%2 [Umax 

n p = ±X. / da ^ (Aii) 



^ = ~^ml (A12) 

and 



f maX da v^arctan ( ^ V ™'f \ (A13) 

^min V ^ Pmax / 



(i) Tp max , . 

?7e = 5— Wm« - W mmJ (A14) 

7T' 3 



( 2) Tp max [*»" 2 2 /lW2$ + $2\ 



V 1 V iv mm T ^mio/ \ max min / 
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and 



Tjp; 



,3 

max 



3V% 

+ 2v / 2($ 2 



V2\n (I 



$4 



+ $ 4 ■ 

1 mm 



2y / 2arctan 



<T>2 -I- <J) 2 . 

max 1 mm 



$2 $2 _ 1 

max mm 



max mm I 



4$ max arGtan 



I v ^ ^ max I 

\ max / 



4$1- arctan 



2 

min 



\ $ 2 • - 1 

\ mm 



(A18) 



We now go over to compute the contribution to r\ c coming from the sound waves. We 
have 



4Tf 
' 3tt 3 



/ Jd\2 



da" (a 



da 



a 



(A19) 



[( (T 0)2_ (T 2/3]2 + (4 7a 2 a 0/3)2 

To provide an analytic result we shall be concerned only with the case of small 7. The idea 
is that the behavior of Eq. (1A19j) for 7 — > is dominated by the singularity at (a ) 2 = a 2 /3. 
To isolate the singular behavior, we write 

4T7 2 



Vc 



3tt 3 



da 



Pmax 



da 



(a°)V 4 



[(a ) 2 - a 2 /3 + 4z 7( T 2 a /3][(a ) 2 - a 2 /3 - 4i7<7 2 <7°/3] 

(A20) 



Completing the squares, the integrand reads 



(a°) 2 a* 



[(a + 2i 7 <x 2 /3) 2 - a 2 /3 + 4 7 V 4 /9] [(a - 2i 7 a 2 /3) 2 - <x 2 /3 + 4 7 2 a 4 /9] 
Since we are interested in the small 7 behavior, we neglect the 7 2 terms so 



(A21) 



4T 7 2 



d<7° 



Pmax 



da 



(a°) 2 a 4 



(A22) 



[(<r° + 2i 7 a 2 /3) 2 - a 2 /3][(a° - 2t 1 a 2 /3) 2 - a 2 /3] 

We assume u m i n is small enough and oj max is large enough that we can compute the integral 
by closing the contour. We choose to close it from above. We thereby pick up the poles at 
2i 7 a 2 /3 ± a/ VS. We get 



4iT 7 2 
3tt 2 



Prna x 



da 



(2z 7 a 2 /3 + a/V3) 2 a 4 



_ [(a/V3 + 4z 7 a 2 /3) 2 - a 2 /3] [2a 2 /3] 
(2i 7 cr 2 /3 - a/V3) 2 a 4 



(A23) 



[(-a/V3 + Ai^/3) 2 - a 2 /3] [-2a 2 /3] . 
We keep only the leading 7 terms, whereby the integral becomes trivial and gives ~ Pmax/l- 



The final result is 



T'jp: 



max 



(A24) 
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Appendix B: Estimating 



In this appendix we describe the procedure used to estimate the value of oo m i n , the 
frequency that divides ultrasoft from soft modes. The idea is to find u min by requiring that 



the soft mode equations can be linearized [56|, 159 ]. 



The condition for a mode with wavevector k to be linearizable is 

k{v 2 ) <uk 2 y/<yj (Bl) 

where 

V = H = ±L (B2 ) 

and (v 2 ) is the thermal average of the square of the soft mode velocity. This thermal average 
can be estimated from the magnitude of the shear correlator S12 given in Eq. (1A6j) . For our 
purpose of estimating bJ m %n it is enough to consider only the shear part of the correlator. We 
get 

= / da o da ° (B3) 

where we have used that s = aT 3 and set co max = p max . For the a integral, we call 

a = ^tan0 (B4) 

so we get 

87r fP max I" f p T \ ( kT \ 1 

^ = ^J k da a2 r cta \wp ) _ arctan l(^pJ J (B5) 

If rj/s is small, the arguments of both arctans are large, and we may approximate 

arctan(-) ps <5 (B6) 



to get 



Putting all together we get the condition for soft modes as 



Pmax J \ S / Pmax 

with y some constant of order one. The boundary lies at 



k \ fv\ 1/2 Tk 

^Vl 1 ) — (B8) 



v * ' Pmax 
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This is our final result. For simplicity, in order to estimate the value of £ we will set y = 1. 
For later convenience, we will rewrite Eq. (IB9I) in terms of % — p ma xV / (sT) . We have 

1 



(BIO) 



' 1 + (2)3/21 

For a strongly coupled theory with x = 1/6 we get from Eq. ( 1B10I) that £ = 0.88 for 
rj/s — 0.08 and £ = 0.72 for rj/s — 0.16. For a weakly coupled theory with \ = 1/12 we get 
£ = 0.78 for rj/s = 0.08 and £ = 0.56 for rj/s = 0.16. 



For completeness, we briefly discuss the approach followed in [52J to estimate the value of 
^min, which is different from the one adopted here. There it was found that the contribution 
of the quadratic term (p + Pq)u^u v to the retarded correlation function of two stress energy 
tensor operators T xy gives rise to a nonanalytic term that goes like co 3 ^ 2 . The value of uj m %n 
can be found by comparing the standard term T]T n u> 2 to the nonanalytic correction. For 
uj < 0J min , second order fluid dynamics ceases to be valid because the gradient expansion 
breaks down, because the nonanalytic correction becomes larger than the standard term. In 



our notation, the result can be expressed as follows (see Eq. 4.8 of [52| ) 



^^r^i (bid 



a 2 \ s 



From Eq. (1B1 lj) we get, for the QGP with a = 10 and x = Pmaxl = 1/6, that £ = 0.42 for 
rj/s = 0.08 and £ = 0.007 for rj/s = 0.16. These values are smaller than the ones we obtain 
from requiring that the soft mode equations be linearizable. 
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